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Abstract 

The conformation and the phase diagram of a membrane protein are inves- 
tigated via grand canonical ensemble approach using a homopolymer model. 
We discuss the nature and pathway of a-helix integration into the membrane 
that results depending upon membrane permeability and polymer adsorptiv- 
ity. For a membrane with a permeability larger than a critical value, the 
integration becomes the second order transition that occurs at the same tem- 
perature as that of the adsorption transition. For a nonadsorbing membrane, 
the integration is of the first order due to the aggregation of a- helices. 

PACS numbers: 87.15.By 64.60. Cn, 61.41. +e 
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Membrane proteins regulate signal transduction and ionic or macromolecular transport 
across biomembranes. Because of their unique roles in biological functions, their confor- 
mations and the folding pathways are important issues in biological physics no less than 
the corresponding aspects of globular protein folding. Folding of integral membrane protein 
carries different characteristics compared to globular protein folding due to the hydrophobic 
environment of phospholipid membrane In watery solvent, the outer surface of glob- 
ular protein is usually covered with hydrophilic segments, while the inner space is filled 
with hydrophobic segments to minimize the protein-solvent interaction energy. In contrast, 
membrane proteins have hydrophobic outer regions inside the membrane to minimize the 
protein- lipid interaction energy 0. 

The three dimensional structures of a great variety of globular proteins are experimentally 
known. Yet the structures of only a few membrane proteins are resolved, because the proteins 
embedded in hydrophobic membrane are difficult to handle experimentally 0. Since the 
structural determination of Bacteriorhodopsin , the idea has been widely accepted that the 
membrane proteins are predominantly made up of a- helices induced by hydrogen bonding . 
Unlike the globular proteins, the membrane proteins can adopt only a few basic structures 
such as a-helix, allowing more tractable theoretical approaches for membrane proteins. 

While a number of theoretical studies have been done separately on globular protein 
folding 1^ and polymer adsorption on membranes 0, there are few efforts devoted to the 
folding of membrane proteins involving the surface adsorption ||^. In this Letter, we address 
this problem using the statistical mechanics via grand canonical ensemble approach. To 
extract the salient features of the conformations and their transitions from the intractable 
complexity characteristic of the real proteins, we introduce a simple but tenable model: 
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a long homopolymer which undertakes a random walk outside the membrane regarded as 
planar, and can interact with it via contact binding on its surface and penetration into its 
interior(Fig. 1), as will be detailed. Motivated by the fact that hydrogen bonding is very 
stable in the hydrophobic environments we assume that a-helix structure is formed if 
and only if the segments are placed within the membrane. Here we neglect other secondary 
structures such as /^-sheets for simplicity. Another important observation to incorporate is 
that the a-helices preferentially aggregate to form a thermodynamically stable structure, 
called a-helix oligomer, which is dominant over the dispersed a-helices 

In our model, an a-helix column has fixed number of hydrogen bonds n (implicitly rep- 
resentative of membrane thickness), with the statistical weight Wh ~ exp(— /3(neft + e^)), 
where P = l/fc^T, at = exp^AsM/ks) < 1- The < and Asm < are the energy and 
the entropy change associated with hydrogen bonding, and Ca < is the aggregation energy 
per helix column. On the membrane surface, polymer segments are allowed to be adsorbed 
with the statistical weight for k segments given by Ws ~ exp{—j3kes), where < is 
the segmental attraction energy, and a = exp{As/kB) < 1 with As the segmental entropy 
change by adsorption. Due to chain connectivity, the domains other than the membrane- 
bound oligomer and surface-adsorbed trains consist of two end tails, loops starting or ending 
at surface trains, and loops that connect two helix columns(return to the starting point as 
the radius of an a-helix will be neglected in this work) (Fig. 1). A tail of k segments has the 
statistical weight of random walk in half space, which departs from the membrane surface 
and ends up anywhere in the half space [rU|, Wt ~ q'^/k^^'^, where q is the segmental par- 



tition function. Loops have the statistical weights of random walks ending up anywhere on 
the surface, Wl ~ q'^/k^^'^, and of those returning to the starting point, Wr ~ q^ /k^^'^ [ITT . 



We construct grand partition functions of five different kinds of domains as follows. The 
bound chains as a helix column and a surface-adsorbed train have respectively the partition 
functions, 

Qh{zm) = IhZM(T^exp{-(3{neh + (1) 

oo 
k=l 

= Iszaexp{-Pes)/{1 - za exp{-Pes)), (2) 

where Ih, Is are the nucleation or initiation parameters for an a-helix and a surface- adsorbed 
train respectively [jl2[. The z and zm are the segmental fugacities outside and inside 



membrane respectively, which defines the chemical potential difference A/x = fiM ~ = 
f3~^\og{zM/z), a measure of membrane permeability determined by environmental effects 
such as membrane hydrophobicity. With q = 1, neglecting the irrelevant bulk contribution, 
the grand partition functions of a tail, a loop, and a returning loop are given by 

oo 

Qriz) = At E ^'/k'^' = ^T9i/2{z) (3) 

k=l 

oo 

QLiz) = At E ^'/k'^' = ALgs/2iz) (4) 

k=l 

oo 

QRiz) = ArY^ z'^/k'/^ = Ang^dz), (5) 
fc=i 

where At, Al, Ar are constants of order of unity, and gm{z) is the polylogarithmic function 
of order m. 

The total grand partition function of the membrane-bound polymer can now be calcu- 
lated considering every possible conformation made of all the domains. To this end, consider 
the transfer matrices defined as 
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Ql 
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which represent respectively the two membrane-bound domains(adsorbed-train and a-hehx) 
and two types of loops(L and R) joining them. Introducing the matrix B = X Y1^=q (YX)^, 
whose elements properly incorporate the interconnected arrays of all the domains between 
two tails (for example, the element B12 represents the connected arrays starting from train 
and ending with helix), we finally get the total partition function Q and its diagrammatic 
representation(Fig. 2), 

Q — Q\'{Bii + B12 + B21 + B22) 
= [QHQs + Qh) 

+ [QHQIQl + QsQlQh + QhQlQs + QIQr) 

+ terms with p > 2. (6) 

Here we regarded the two ends of the polymer to be distinguishable. We will consider the 
thermodynamic limit, written as {N) — > 00, in order to define phase transitions that result. 
Then the partition function is reduced to 

Qti -^c 

(i-A)-i, T<re, 

where A is the largest eigenvalue of YX, given by 

1/2 



(7) 



A = \{.QsQl + QhQr) + \ [{QsQl - QhQr? - ^QsQhQI 



(8) 



and Tc is a critical temperature determined from X{z = 1, T = Tc) = 1. 

For T > Tc, tail is the only allowed conformation, indicating that the polymer tends to 
be desorbed to the outer region of the membrane. For T < Tc, the surface trains and/or 

5 



helices with loops and returning loops in between becomes dominant, which indicates the 
stability of membrane-bound phase. The segmental fraction of each domain can be de- 
fined as fi = {Ni)/{N) = {Ny\d\ogQ/d\ogQi){d\ogQi/d\ogz), where i = S,R,L 
represent three different types of domains considered(Fig. 1), and fn = {Nh) / (N) = 
{N)~^ {dlog Q/dlogQH){dlogQH/dlog zm), for helix domain. Depicted in Fig. 3 are the 
segmental fractions vs. temperature, where all the energy parameters, scaled in units of 
le^l, are taken to be the same order of magnitude, and the entropies As and Asm are taken 



to be order of ks |T3[- As shown in Fig. 3(a), the desorpt ion- adsorption transition, which 
is of the second order as is known, takes place at T = Tc, where the order parameter, the 
surface-adsorbed segmental fraction(/5'), increases from zero continuously as T is lowered 
from Tc. Further lowering of the temperature drives the polymer integration in a form of 
a-helix aggregate at T = Th, the helix inclusion temperature, which is defined by the local 
maximum of specific heat. The specific heat curve, shown in Fig. 4, clearly indicates that 
structural changes occur at both Tc and Th- The specific heat diverges at T = Tc, and has 



a local maximum at T = T/^ indicating the helix formation that is a crossover ||14|| . 

As the value of = e/^ — A/i, the energy of a-helix inclusion per segment, is lowered, or 
membrane permeability is increased, Th approaches Tc, so that for the values of smaller 
than a critical value (about —2.6 using the parameters employed in Fig. 3), the helix inclusion 
is promoted to the second order transition with Th and Tc. Figure 3(b) depicts the segmental 
fractions for this case(e^ = —3). It is shown that, in contrast to Fig. 3(a), helix formation 
dominates over adsorption below the common transition temperature. The phase diagram 
in Fig. 5 summarizes the foregoing discussions concerning the conformational phases and 
their transitions for wide range of ejj and temperature. 



The two temperatures appearing in our model, the desorption-adsorption temperature 
Tc and the hehx inclusion temperature T/j, govern the pathway of our model- homopolymer 
integration in membrane. They are respectively similar to coil-to-globule transition temper- 
ature Te and folding transition temperature Tf in globular protein folding; like the globular 
phase, the adsorbed phase is indeed an intermediate state approaching the native folded 
structure Recently Klimov and Thirumalai showed an evidence that globular protein 



folding time is scaled as r ~ exp [J\Tg — Tf\/Tg], where J is a model-dependent constant 
T6| . Even without considering the analogy, the Fig. 5 suggests that a transition into rapid 



a-helix integration in membrane can be attained for permeable and adsorbing membranes, 
with lower than a critical value where = T^. A detailed analysis of the free energy 
landscape and barrier crossing dynamics should confirm this highly plausible suggestion. 

We now consider the situation of nonadsorbing membrane where Qs-, the partition func- 
tion of surface-adsorbed domain vanishes (for example, /<,, the adsorption initiation param- 
eter is zero). We then find the eigenvalue A = QhQr, signifying that the partition function 
incorporates the conformations generated from recurrence of a helix and a returning loop 
in series. In this case, the polymer inclusion forming the a-helix oligomer is found to be 
the first order transition. At the transition temperature T^, which is again determined by 
\{z = 1, Tfi) = 1, the order parameter, the fraction of segments in the a-helix oligomer (/j^), 
changes discontinuously from zero to 

AfH = n/[n + gy2{l)/g,/2{l)], (9) 

which incurs the segmental latent heat of dissociating the membrane-bound oligomer, 
Ih = l^l + ea/n\AfH. The reason why the transition should be discontinuous in the ab- 
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sence of adsorption is argued as follows. The inclusion accompanies the aggregation of 
helices, which restricts the loops to be closed between helix columns. Presence of long re- 
turning loops is suppressed entropically, and, also energetically in favor of the aggregate, 
which tends to be of significant fraction at the transition from the desorbed phase. This is 
in sharp contrast to the two second order transitions that can be obtained from our theory, 
the desorption-adsorption transition without helix inclusion(A = QsQl) and the inclusion 
transition in forms of dispersed helices without aggregation and adsorption(A = QhQl with 
= 0), where the loops can be transformed continuously into the adsorbed segments and 
helix columns respectively. Similar discontinuous transition was reported in the helix-coil 



transition of double strand DNA [11,17]. For the biological processes without appreciable 



changes in temperature and volume, certain mechanism of latent heat involving enzymatic 



activity could be essential to facilitate the first order transition |Tl]]. Furthermore, the con- 
clusions of this and foregoing paragraphs, if properly extended to an asymmetric membrane 
where one side is adsorbing and the other is not , imply that the rapid protein integration 
can be promoted only through the adsorbing side. It would be important to confirm this 
possibility as well as our results for symmetric membranes by experiments and/or simula- 
tions. 

In summary, we studied various membrane-protein conformations and different pathways 
to the native structure of an a-helix aggregate as a function of temperature and membrane 
permeability. Two significant conclusions obtained are: 1) The polymer inclusion pathway 
is determined by membrane permeability above a critical value of which the adsorption and 
the helix inclusion converge and 2) The nature of inclusion transition is determined by the 
availability of the polymer adsorption on membrane surface, due to the chain connectivity 
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constraint. Although the important sequential heterogeneity and finite length effect involv- 
ing chain stiffness of real proteins are neglected, our model gives some nonspecific features 
of membrane-protein conformation, in particular, the roles of membrane hydrophobicity and 
segmental interaction with the membrane surface. 

We acknowledge the support from KOSEF(961-0202-007-2), BSRI(97-2438), and 
POSTECH special fund program. 
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FIGURES 



FIG. 1. Schematic figure of a membrane protein. Five different domains are indicated. 

FIG. 2. Diagrammatic representation of the partition function. 

FIG. 3. Segmental fractions of a membrane-bound polymer vs. temperature. The parameter 
values throughout this paper are selected as[13] a = exp(— 1), = exp(— 2), Al = 1, Ar = 0.5, 
Ifi = 0.01, Is = 0.1, n = 12, £„ = —3, where energies and temperature are in units of |es|. (a) 
e\ = Eh — A|Lt = —2. Desorption-adsorption transition takes place at T = Tc ~ 1.30 continuously, 
while the helix inclusion occurs at T < T/j ~ 0.82. (b) ejlj = — 3(more permeable membrane). 
Adsorption and helix inclusion occurs simultaneously at T = Tc = ~ 1.46, but helix structure 
becomes dominant over adsorbed state at lower temperatures. 

FIG. 4. Specific heat (in units of ks) versus temperature. Parameters are the same as in 
Fig. 3(a). 

FIG. 5. Phase diagram of a membrane-bound polymer. Solid line indicates the second order 
transition, while dotted line the crossover between adsorption and a-helix inclusion. Parameters 
are the same as in Fig. 3. 
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